BCS Theory for Binary Systems with 2D Electrons 
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MgB2 is considered as a binary system with 2D electrons. The classic BCS theory is applied to this system. 
The transition temperature (Tc) is found to be relatively high, because 2D electrons are more capable of moving 
with the atoms, on top of other special features of this system to enhance the electron-phonon interaction. This 
system may also shed light on the nature of superconductivity in cuprates. 



PACS numbers: 74.20. Fg, 74.25.Bt, 74.80.-g 

There is little doubt that MgB2 (T^ = 39K) is a conven- 
tional superconductor ||l|. The isotope effect is found to 
be strong, and the energy gap is also of the conventional 
type 0, pi. Unsolved is the mystery about Tc which 
seems to be either above or at the limit suggested by 
the theory of Bardeen, Cooper and Schricffer (BCS) for 
phonon-mediated superconductivity [Q. Therefore it is 
of current interest to apply the BCS theory to MgB2, 
taking the special features of this compound into con- 
sideration, to see if a proper Tc would arise. 

MgB2 features alternating B and Mg rings ||l[. En- 
ergy band calculation reveals that superconductivity in 
MgB2 is driven by the 2D bounding a bands |^. We 
consider a simple model of alternating layers of light 
and heavy atoms, to which electrons pertain (Fig, 1). 
These electrons are more capable of moving with the 
atoms, compared with free electrons, to enhance the 
electron-phonon interactions. These interactions may 
take place either within the same atom layer (intra-layer 
coupling), or between neighboring layers (inter-layer cou- 
pling). The intra and inter-layer couplings may join force 
to raise Tc still further. 

Our simple model might help us to understand cup- 
rates, regardless of the nature of the pairing mechanism. 
In cuprates atoms are also layered, with electrons per- 
taining to the atom layers [p[. The intra and inter-layer 
couplings are in competition: if none dominate then both 
are suppressed. This leads to a forbidden zone of dop- 
ing, where the inter-layer coupling is stronger than the 
coupling within the heavy atom layers, but weaker than 
that within the light atom layers, so that superconduc- 
tivity vanishes, reminiscent of the observation by Mood- 
enbaugh et al. in La2-i;Baa;Cu04 Q]. Careful doping 
may leave singlet pairs inside the forbidden zone, triplet 
pairs outside: p-wave pairing is possible in the BCS the- 
ory, and indeed in a large number of candidate theories, 
where p-wave symmetry is not intrinsically favored H . 

We apply the BCS theory to our binary system. The 
derivation is simple in principle but involved technically. 
Second quantization is actually based on the assumption 
that, despite its statistical nature, the displacement of a 
particle in a crystal is the exact replica of the displace- 
ment of its neighbors, save a phase factor, which serves 



as a parameter to categorize particle displacement into 
groups known as Fourier components. The operation 
to add or take away a component is symbolized by the 
creation or destruction operator. In this respect there 
is little difference between normal metals and our bi- 
nary system. Phonons (atom displacements) have long 
been treated in complex systems ||l^ . We treat electrons 
with the understanding that proper wave functions can 
be found to reflect the physics of 2D electrons: they per- 
tain to the atom layers, and their wave functions overlap 
with those in neighboring layers. When T > we find 
the following ensemble average of electron energy 
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Fig. 1 Crystal of light (wt. Mi in layer 1, open circles) 
and heavy (wt. A/2 in layer 2) atoms, a, b and c are 
lattice constants, M1/M2 — 0.7, a = b = 0.5c. The anti- 
nodes of ip^ and 4'k s-re in layer 1 and 2, respectively. 
Values of V are in an arbitrary unit. 



Here k is the wave vector to specify the electron state, 
Ck state energy relative to the Fermi level, ft.k occupation 
probability of state k, f^' excitation probability of the 
ground pair at states k and — k, and i and j indicate 
atom layers. Following BCS M we associate k and — k 
with spin | and J, respectively unless stated otherwise. 
We have 
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where uji is the phonon frequency, k = k — k' phonon 
vector, I identifies phonon branches. We have 2 atoms 
in the primitive cell, giving 6 branches in each polariza- 
tion. However, transverse phonons do not interact with 
electrons 10 so that Eq. (g) involves just 6 branches (3 



(«j) 



of them optic). The matrix element A^^ measures the 
strength of electron-phonon interaction, which couples 
electrons as long as their wave functions overlap: the 
coupling may take place either within the same atom 
layer (i — j) or between neighboring layers (i ^ j). We 



have 7W|^'^^ = Mf-^\ so that V^ 
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In the binary system the entropy of the electron en- 
semble is 
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where the factor 4 indicates the degree of degeneracy: at 
the same energy level there can be 2 electrons in each of 
the 2 atom layers. We have assumed the same density 
of states in both atom layers, otherwise the factor 1/2 in 
Eq. (P) must be adjusted. We minimize the free energy 
F ^W -TS with respect to h'^^ and find 
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which are counterparts of the BCS self-consistent equa- 
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F with respect to f^'' and find 
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We also minimize 
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where z = 1 or 2, i.e. Ai = A2 = A, which leads through 
Eqs. (|), (|) and the famihar BCS algorithm [g to 
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k runs over a range where ^huj < Ck < 'hu, lu being 
the average phonon frequency, Vij = {V^^{^, )av matrix 
element averaged over the above range of k. We find 
from Eqs. (0) and (|) that 
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We have assumed M2 > Mi (Fig. 1) so that V22 < Vn 
(to be justified later). When ^22 < V"i2 < (1^11^22)^/^, 
the right hand side of Eq. (ph becomes negative, which 
means /^ > 1/2 must hold at least for some states. 

This not only violates the nature of f^ but also leads 
to a contradiction: A may become so large that Eqs.^) 
and (P) cannot balance. Similarly, we find from Eq. dlC ) 
that (ViiV22)^/^ < V12 < Vli is not allowed. Apparently 



V22 < V12 < Vn 
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is a forbidden zone: if the inter-layer coupling is stronger 
than the coupling within the heavy atom layers, but 
weaker than that within the light atom layers, then su- 
perconductivity is suppressed. Outside this forbidden 
zone we can add Eqs. (^) and (10) together and find 
through Eq. (0) that 
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where the summation over k has been converted into 
integration, iV(0) is the density of states at the Fermi 
surface. Eq. (113) is formally identical to the familiar 
BCS equation W with 
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which is the counterpart of the BCS average matrix ele- 
ment g. 

Now we evaluate V^ j.' ™ ^'i- (0)- ^^1 ^^^^ point we 
have to be more specific about the electron wave func- 
tions. For simplicity, we consider super-positions of two 
Bloch functions: 
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where k — xfej, -I- yky is a 2D wavevector, a spin, c 
lattice constant in the c-direction. Both ip'^^^ and ^^^"^ 



have the same Bloch energy, and a =t or J,, so that 
our electron wave functions are of quadruple degeneracy, 
consistent with the factor 4 in Eq. (||). Furthermore "0*^^^ 
and V'^^'' are orthogonal to each other: they are proper 
base functions for second quantization. 

Both ^l}^^' and '(/'"' are traveling waves in the a-h 
plane but standing waves in the c-direction: they are 
mobile only in 2D. Apparently V' pertains to layer 1 
(its anti-nodes are in that plane) whereas ip^"^' pertains to 
layer 2 (Fig. 1). It is implied that electron waves overlap 
with those in the neighboring layer, but not with those in 
the next neighboring layer (there is a node in between). 
It appears that ^'^^^ and ?/'^^^ are the leading and perhaps 
dominating terms of the electron configuration in a real 
binary system: dispersions arising from Eqs. ( |l4|) and 
(|l^) are similar to that in Q . 

We average V^^, according to the text below Eq. (||) 
and find 
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V22 also is given by Eq. (|l^) with Mi and M2 inter- 
changed. Apparently Vu > 1^22 when Mi < M2. Here 
we have used the Debye phonon model for simplicity, A 
and B are values of the overlap integral function F{x) — 
3(sina;- a;cosx)/a;3, with x = 3Ma^/^Z'^/^C{l - C)^/^ 
and 3.84ai/3[Z2/3C'(l - (') + 0.65(a/c)4/3]i/2^ respec- 
tively, C = 3Zma^{6V)yN{0)niveF, Z is the average 
valency, ep the Fermi energy. Eqs. (|6|) and ( p7| ) exist in 
the sense of the Cauchy principal value (used by Kuper 
to verify the ECS theory) Q. 

The original BCS theory assumes a spherical Fermi 
surface, which may not be the case in binary systems 
P, m. On the other hand, the states of 2D electrons 
form merely a cross-section of the Fermi sea: so long 
as this cross-section is circular, it makes little difference 
whether the Fermi sea is spherical or cylindrical, pro- 
vided that N{0) is of a proper local value. Indeed, all 
the 2D calculations over this circular cross-section par- 
allel the calculations over the spherical Fermi sea: no 
assumptions by BCS have been violated. 

We let M1/M2 = 0.7 and a = 6 = c/2, and evaluate 
Eqs. (|lq ) and (17) numerically (see Fig. 1). The forbid- 
den zone is bordered hy Z = 0.1311 and 0.1359: Vu, 
V12 and V22 are comparable when the valency is rela- 
tively small. This arises from Eq. (y), where the matrix 

(12). 
element Ai, ' is proportional to the component wave 



number n/c in Eqs. dl4) and (Raj. In contrast ^A 



(1,1) 



(or A^ 
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are proportional to kx and ky, which become 



smaller the smaller the Fermi surface {Z smaller). These 
lead to the following exotic properties: 

Critical temperature: which can be relatively high 



for two reasons. First, the intra and inter-layer cou- 
plings may join force to enlarge the energy gap. For 
example, when V^i2 ~ Vu in Eq. (|l3|), V ~ 2Vii holds 
in our binary system, compared with a normal metal 
or alloy where V = Vu. This is hardly a surprise, be- 
cause in Eq. (H^) the interaction Hamiltonian always low- 

(i 7) 

ers the energy of the system, as long as V^^ ^, > and 

/k "^ -'^Z^' regardless of whether the electron-phonon 
interaction takes place within the same atom layer or 
between neighboring layers. Second, electrons pertain- 
ing to atom layers are more capable of moving with the 
atoms. This can be seen from Eq. (16) where Vu is 
proportional to 2A'^/AIi when A = B. In a binary al- 
loy with free electrons Eq. ( p^ also applies, with both 
A + B and A — B being replaced by A. Therefore Vn 
becomes proportional to A^/Mi when Mi = M2: free 
electrons are less capable of moving with the atoms, so 
that the electron-phonon interaction is weakened by half. 
Indeed, in Fig. 1 the value of V for 2D electrons is 3.35 
times larger than that for free electrons when F12 — Vn. 
We use the Debye phonon model when deriving Eqs. dlfl) 
and ([l7|), so that the contribution by optic phonons is 
somewhat overestimated. However, values of V for free 
electrons are also overestimated, so that we can still use 
the data in Fig. 1 to compare the energy gap in a normal 
metal or alloy with that in our binary system. 

Anomalous suppression of superconductivity: which 
is considered as one of the long-standing mysteries asso- 
ciated with high Tc cuprates |12|. It was found that in 
La2-xBa:rCu04 T^ « 25K when x = 0.09 and x = 0.15, 
whereas Tc ~ 5K between these two maxima |7| . Accord- 
ing to energy band calculation, in LaBaCu04 the Fermi 
surface is crossed by the CuO band and two Oz bands 
||6|] : superconducting carriers are in atom layers which on 
average have two different masses. Perhaps LaBaCu04 
could be modeled as a binary system, which falls into 
the forbidden zone when the valancy is adjusted into a 
proper value through doping. 

Spin triplet pairs: which may arise when spin sin- 
glet pairs fall into the forbidden zone. Specifically, while 
V12 < Vn holds for singlet pairs (inside the forbidden 
zone, suppressed), F12 > ^11 may hold for triplet pairs 
(allowed), because Vn declines faster than V12 when the 
pair symmetry changes, due to the stronger effect of the 
exchange term on Vn, which arises from intra-layer cou- 
plings, where electron waves overlap to a greater extent, 
compared with the inter-layer coupling. It is interesting 
that magnetic excitations are observed in YBa2Cu308 
and Bi2Sr2CaCu208 [|3|- According to Yin et al. these 
excitations may not be attributed to the normal state 
properties of cuprates, because the excitation intensity 
decreases continuously with increasing temperature, and 
vanishes above Tc |13| , |l4| . Furthermore, magnetic ex- 
citations are observable only when doping is optimized 
|13|. Perhaps these cuprates may also be modeled as a 



binary system, where careful doping may drive singlet 
pairs into the forbidden zone, leaving triplet pairs out. 
It is difficult to dope MgB2, so that singlet pairs cannot 
be discriminated: the energy gap must be of the BCS 

type §. 

In conclusion, the BCS theory may explain the rela- 
tively high transition temperature of MgB2 which, ac- 
cording to Bud'ko et al., nearly double the previous 
record for a nonoxide and non-C60-based compound [^ . 
We find that in a binary system with 2D electrons V can 
increase 2-3 times (Fig. 1). According to Yin et al. many 
experimental observations in cuprates are commonly fit 
with a phenomenological model that derives from the 
original BCS theory, characterized by a gap equation 
corresponding to a presumed symmetry of the order pa- 
rameter, with an adjustable dimcnsionless coupling con- 
stant, and a given Fermi surface ||lj| . We find that in our 
binary system the p-wave symmetry may arise in a nat- 
ural manner from the BCS theory. This may help us to 
develop a microscopic theory for cuprates. In this letter 
we have used the original BCS theory, where damping 
and retardation are ignored with negligible effect [|l5J. 
This effect is likely to be similar, whether the supercon- 
ductor is a normal metal, or a binary system with 2D 
electrons. Therefore our major conclusion, i.e. V can be 
2-3 times larger in the later case, is likely to stand. In 
future we would like to improve our results by taking 
damping and retardation into consideration. 
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